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We study a system of nonintersecting (1 + l)-dimensionaI fluctuating elastic interfaces ('vicious 
bridges') at thermal equilibrium, each subject to periodic boundary condition in the longitudinal 
direction and in presence of a substrate that induces an external confining potential for each interface. 
We show that, in the limit of a large system and with an appropriate choice of the external confining 
' potential, the joint distribution of the heights of the A'^ nonintersecting interfaces at a fixed point on 

^— »' , the substrate can be mapped to the joint distribution of the eigenvalues of a Wishart matrix of size 

A'^ with complex entries (Dyson index f3 — 2), thus providing a physical realization of the Wishart 
' matrix. Exploiting this analogy to random matrix, we calculate analytically (i) the average density 

^ , of states of the interfaces (ii) the height distribution of the uppermost and lowermost interfaces 

(extrema) and (iii) the asymptotic (large A'^) distribution of the center of mass of the interfaces. In 
' the last case, we show that the probability density of the center of mass has an essential singularity 

around its peak which is shown to be a direct consequence of a phase transition in an associated 
Coulomb gas problem. 

1^ ' PACS numbers: 05.40.-a, 02.50.-r, 05.70.Np 
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^ . The system of N nonintersecting elastic lines was first studied by de Gennes [l| as a simple model of a fibrous 
^ ' structure made of (1 + l)-dimensional nonintersecting flexible chains in thermal equilibrium, under a unidirectional 
I stretching force. These elastic lines can also be viewed as the trajectories in time of N nonintersecting Brownian mo- 

ations, a system studied in great detail by Fisher and co-workers context of commensurate-incommensurate 

I ^ (C-IC) phase transitions. In this context the nonintersecting lines are the domain walls between different commen- 
i surate surface phases adsorbed on a crystalline substrate. The 'nonintersection' constraint led Fisher to call this a 
' problem of 'vicious' random walkers who do not meet (or kill each other when they meet). Since then, the vicious 
O ,. walkers model has had many physical applications, e.g., in wetting and melting 0, [sj, as a simple model of polymer 

1 ^ I ■ network Q , in the structure of vicinal surfaces of crystals consisting of terraces divided by steps Q and also in the 

[ context of stochastic growth models Q . 
04 ' Depending on the underlying physical system being modelled by N vicious walkers, one can pose and study a 
^ ' variety of statistical questions. For example, Huse and Fisher [2] studied the so called 'reunion' probability, i.e., 
\l . the probability P(t) that N vicious walkers starting at the same position in space reunite exactly after time t at 
their same initial position but without crossing each other in the time interval [0, t] and showed that it decays as a 
power law P{t) ~ for large t. Another pertinent issue is: given that the N walkers have reunited for the first 

' time at time t, what can one say about the statistics of the transverse fluctuations of the positions of the walkers 
i at any intermediate time < r < t (see Fig. 1)? Such configurations where N walkers emerge from a fixed point 
in space and reunite at the same point after a fixed time t are called 'watermelons' as their structure resembles 
that of a watermelon (see Fig. 1). Such a watermelon configuration also describes the structure of the 'droplet' or 
^ ' the elementary topological excitation (vortex-antivortex pair) on the commensurate (ordered) side of a C-IC phase 
. ^ transition P, (sj with the longitudinal distance between the pairs being t. The statistical (thermal) fluctuations of the 
transverse sizes of such watermelon defects play an important role near the phase transition. This initiated a study 
^ ] of the tranverse fluctuations of the nonintersecting lines in the watermelon geometry with fixed longitudinal distance 
t. In the random walk/probability language, this means studying the transverse fluctuations of the trajectories of N 
walkers conditioned on the fact that they started and reunited at the same point in space after a fixed time t without 
crossing each other in between. Another similar interesting geometrical configuration is a 'watermelon with a wall', 
i.e, N nonintersecting walkers starting and reuniting at the same point in space (say the origin) after a fixed time t, 
but staying positive in [0,t] (see Fig. 1). 

Recently, the transverse fluctuations of N nonintersecting lines have been studied extensively in watermelon ge- 
ometry over [0, t\ both with and without a wall and important connections to random matrix theory have been 
discovered 0; M S [13, [Hi IHj 13 ■ For example, the joint distribution of the positions of all the walkers at a fixed 
time < T < t for watermelons without a wall was shown to be identical (after appropriate rescaling) to the joint 
distribution of eigenvalues of a Gaussian random matrix belonging to the unitary (CUE) ensemble [a. [ll|. On the 
other hand, the joint distribution of the positions at a fixed time < r < t for watermelons with a hard wall at the 
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FIG. 1: Typical watermelon configurations without a wall (left) and with a wall (right) at the origin for N = 5 nonintersecting 
Brownian motions over the time interval < t < t. 



origin was computed recently |13| and was shown to be identical (after an appropriate change of variable) to the joint 
distribution of eigenvalues of a random matrix drawn from the Wishart (or Laguerre) ensemble at a special value of 
its parameters, which also corresponds to the chiral Gaussian unitary ensemble of random matrices [lOj. 

It is useful at this point to recollect the definition of a Wishart matrix. A Wishart matrix W is an {N x N) square 
matrix of the product form W = X where X is a {M x N) rectangular matrix with real or complex entries and 
X'^ is its Hermitian conjugate. If the entries Xij represent some data, e.g., Xij may indicate the price of the j-th 
commodity on the i-th day, then W is just the (unnormalized) covariance matrix that provides informations about the 

correlations between prices of different commodities. If X is a Gaussian random matrix, P{X) cx exp —^Tt{X^X) 

where the Dyson index (3 = 1,2 corresponds respectively to real and complex matrices, then the random covariance 
matrix W belongs to the Wishart ensemble named after Wishart who introduced them in the context of multivariate 
statistical data analysis [l5l |. Since then the Wishart matrix has found numerous applications. Wishart matrices 
play an important role in data compression techniques such as the "Principal Components Analysis" (PCA). PCA 
applications include image processing [l^ [l^, [il|jbiological microarrays [l^, US] > population genetics [U , 
finance 0, H^, meteorology and oceanography The spectral properties of the Wishart matrices have been 

studied extensively and it is known that for M > N, all N positive eigenvalues of W are distributed via the joint 
probability density function (pdf) 

P^.(Al, A^) = e-^ ^ J] |A, - A,f (1) 

k=l i<j 
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where Km is a normalization constant and the Dyson index /3 = 1,2 (respectively for real and complex X). In the 
"Anti-Wishart" case, that is when M < N , W has Al positive eigenvalues ( and N — M eigenvalues that are exactly 
zero ) and their joint probability distribution is simply obtained by exchanging M and N in the formula ([1]). Hence 
we will focus only on the Wishart case with M > N. Note that even though the Wishart pdf in Eq. ^ was obtained 
for integer M > N, the pdf is actually a valid measure for any real continuous M > N . In particular, for Ad — N = ^ 
(M — N — — i) in the case P — 2, the pdf in Eq. ([1]) is realized as the distribution of the squares of positive 
eigenvalues of class C {D) random matrices |lOj, in the classification of Altland and Zirnbauer '^2§\. In addition, the 
pdf in Eq. ([T]) is also realized for d-dimensional squared Bessel processes under nonintersection constraint with [3 = 2 
and M - N ^ d/2 - 1 for d > [2^. 

In Ref. 13], the joint pdf of the positions of nonintersecting Brownian motions in the 'watermelon with a wall' 
geometry, mentioned in the previous paragraph, was shown to correspond to the joint pdf of Wishart ensemble in 
Eq. Il]) with special values of the parameter (3 — 2 and M — N — 1/2. At these special values, the joint pdf also 
correspond to those of the squares of the eigenvalues of class C matrices Jii] . A question thus naturally arises whether 
it is possible to find a nonintersecting Brownian motion model that will generate a Wishart ensemble in Eq. ([1]) with 
arbitrary values of the two parameters /3 and M — N >Q. In this paper we address precisely this issue and show how 
to generate the Wishart ensemble with arbitrary positive (3 and M — N > 0, starting from an underlying microscopic 
model of nonintersecting Brownian motions. 

In this paper we study the nonintersecting Brownian motions in a geometry different from that of the watermelons 
discussed above. Here wc consider a system of (1 + l)-dimensional nonintersecting fluctuating elastic interfaces 
with heights hi{x) [i = 1,2,3, ...A) that run across the interval x E [0, L] in the longitudinal direction (see Fig. 
2). Equivalently the heights {hi{x)} can be thought of as the positions of A nonintersecting walkers at 'time' x. In 
contrast to the watermelon geometry, the lines here are not constrained to reunite at the two end points. Instead, each 
line satisfies the periodic boundary condition in the longitudinal direction, i.e., they are wrapped around a cylinder of 
perimeter L (see Fig. 2). In addition, there is a hard wall (or substrate) at h = that induces an external confining 
potential V{hi) on the i-th interface for all 1 < i < A'. In a slightly more general version of the model, one can also 
introduce a pairwise repulsive interaction between lines. In presence of the external confining potential V{hi), the 
system of elastic lines reaches a thermal equilibrium and our main goal is to compute the statistical properties of the 
heights of these lines at thermal equilibrium. More precisely, we compute the joint distribution of heights of the lines 
at a fixed position < cc < L and show that, in the limit L — > oo (large system), this joint distribution, after an 
appropriate change of variables, is precisely the same as the joint distribution in the Wishart ensemble. Note that 
due to the translational symmetry in the longitudinal direction (imposed by the periodic boundary condition), this 
joint pdf of the heights is actually independent of x. 

Thus our model is actually closer to the solid-on-solid (SOS) models at thermal equilibrium in presence of a 
substrate [30j . except with the difference that here we have multiple nonintersecting interfaces. This model is thus 
appropriate to describe the interfaces between different co-existing 'wet' phases of a multiphase two-dimensional fluid 
system on a solid substrate or a film We note that nonintersecting Brownian motions in an external harmonic 
potential was studied recently by Bray and Winkler [3l| , but they were mostly interested in calculating the probability 
that such walkers all survive up to some time t. In the Brownian motion language, we are here interested in a different 
question: given that each walker survives up to 'time' L and comes back to its starting position, what is the joint 
distribution of the positions of the walkers (or equivalently the heights of the interfaces) at any intermediate 'time' 
0<x <L7 

In this paper, we consider the external confining potential of the form 

V{h) = + ,^ ' with & > and a > 1 (2) 

2 2n,2 

with a harmonic confining part and a repulsive inverse square interaction. Such a choice is dictated by the following 
observations. The harmonic potential is needed to confine the interfaces as otherwise there will be a zero mode. The 
repulsive inverse square potential has an entropic origin. For a single interface near a hard wall, Fisher Q indeed 
showed that the effective free energy at temperature T behaves as ksT/h'^ where h is the distance of the interface (or 
the walker) from the wall. Thus it is natural to choose the external potential of the form as in Eq. In addition, 
as we will see later, such a physical choice also has the advantage that it is exactly soluble. We will see indeed that 
this choice of the potential generates, for the joint density of heights at a fixed point x and in the limit of a large 
system (L oo), a Wishart pdf in Eq. ([1]) with fixed (3 = 2, but with a tunable AI — N = a — 1/2 where a sets the 
amplitude of the repulsive inverse square potential in Eq. Thus when a — 1/2 = A; is a positive integer fc, this 
generates a physical realization of Wishart matrices with integer dimensions AI = N + k and A^. 

Our model is also rather close to the realistic experimental system of fluctuating step edges on vicinal surfaces of 
a crystal in presence of a substrate (or hard wall). When a crystal is cut by a plane which is oriented at a small 
nonzero angle to the high-symmetry axis, one sees a sequence of terraces oriented in the high-symmetry direction that 
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FIG. 2: Nonintersecting fluctuating interfaces with heights hi{x) for < x < L with periodic boundary conditions, hi{0) = 

h, {L). 

are separated by step edges which can be modelled as 'elastic' nonintersecting lines or trajectories of nonintersecting 
Brownian motions [5|. In an external confining harmonic potential but in absence of a wall at h — (such that 
h — > —h symmetry is preserved), the joint distribution of the heights of the lines at equilibrium can be mapped to the 
GUE ensemble @ , although most studies in this context are concerned with the so called Terrace- Width distribution, 

i. e., the distribution of the spacings between the lines. In our model, due to the presence of the wall which breaks 
the h —h symmetry, new interesting questions emerge. For example, it is natural also to ask for the distribution 
of the minimal (maximal) height, i.e., the height of the line closest (farthest) from the wall. Using the mapping to 
the Wishart random matrix, the minimal and maximal height correspond respectively to the smallest and the largest 
eigenvalue of the Wishart random matrix. In addition, it is also interesting and physically relevant to investigate 
the statistics of the center of mass of the N nonintersecting Brownian motions. We will see that strong correlations 
between the lines violate the central limit theorem resulting in strong non-Gaussian tails in the distribution of the 
center of mass. 

Let us summarize below our main results. 

• Using a path integral formalism we show that the computation of the equilibrium joint distribution of heights at a 
fixed point in space can be mapped to determining the spectral properties of a quantum Hamiltonian. Subsequently, 
when the external potential is of the form in Eq. this quantum potential turns out to be integrable and allows us 
to compute the equilibrium joint distribution of heights exactly. In particular in the limit of a large system (L — s- oo), 
we show that the joint distribution of heights, under a change of variables = Ai, is exactly of the Wishart form 
in Eq. ^ with parameters (3 = 2 and M — N = a — 1/2 where a appears in the potential in Eq. ([2|). Knowing the 
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exact joint distribution, we then compute various statistical properties of the heights of the interfaces as hsted below. 

• We find that the average density of lines at height /i, in the limit of a large number N of interfaces, is a quarter 

of ellipse as a function of h, with finite support over [0, 2y^] where b appears in Eq. ([2]). The typical height thus 

scales with TV for large TV as ft-typ ~ VlV- This differs considerably from the case of TV interfaces that are allowed to 
cross, where the typical height is of order one. The spreading of nonintersecting interfaces is a consequence of the 
strong interaction between them induced by their fermionic repulsion. 

• We study the height distribution of the topmost (farthest from the wall) interface in the large TV limit. We show 
that the average height of the farthest interface (maximal height) is 2-^/^ for large TV : it is given by the upper bound 

of the average density of states. The typical fluctuations of the maximal height around its mean are distributed via 
the Tracy- Widom distribution [H, [H, [s^l • However, for finite but large TV, the tails of the distribution of the maximal 
height show significant deviations from the Tracy- Widom behavior. We compute exactly these large deviation tails. 

• We also study the statistics of the height of the lowest (closest to the wall) interface (minimal height) and argue 
that, for large TV, it scales as TV~^/^. This should be compared to the case of non-interacting Brownian motions 
where the typical distance of the closest (to the wall) walker is of ~ 0(1) from the wall. This is again an effect of the 
strong interaction between the interfaces: their mutual repulsion pushes the lowest interface closer to the substrate. 
We further show that the full distribution of the minimal height can be exactly computed for a special value of the 
parameter a = 3/2 in Eq. ([2]). 

• Finally we study the distribution of the center of mass of the heights Gn = ^1+-+''" for large TV. Thanks to 
the analogy between the Wishart eigenvalues and a Coulomb gas of charges, the mean and variance of the center of 
mass can be computed, as well as the shape of the probability distribution: we show that the pdf of Gn, P{Gn — v), 
has a non-analytic behavior (essential singularity) at v = (Gn) (which is shown to be a direct consequence of a phase 
transition of 'inifinite' order in the associated Coulomb gas problem). In addition, we find exact asymptotic results, 

to leading order for large TV, for the mean (Gn) = i^^y^ and the variance (G^) — (Gjv)^ ~ T^nB- 

The rest of the paper is organized as follows. In section|lTl we present our model, compute (via path integral method) 
the joint probability distribution of the heights of the interfaces and compare it to the probability distribution of the 
eigenvalues of a Wishart matrix. In section Hill we analyse some statistical properties of the model. We first present 
the results for the average density of lines (subsection IIII A|) . We then compute the behavior of the maximal height 
(jlll B|) and the minimal height (jlll Cp . Finally we study the distribution of the center of mass of the heights in (jlll Pp . 
Section ITVl concludes the paper with a summary and outlook. 

II. THE MODEL 

Our model consists of TV nonintersecting (1 -t- l)-dimensional interfaces over a substrate of size L (that induces 
an external potential). For simplicity, we first present the model for a single interface in subsection III Al and show 
how, using a path-integral formalism, one can map the problem of calculating the equilibrium height distribution of 
the interface to computing the spectral properties of a quantum Hamiltonian. In particular, calculating the height 
distribution in the limit L — > 00 (large system) corresponds to calculating the ground state wavefunction of this 
quantum Hamiltonian. Then we present the interacting model for general TV interfaces in subsection IIIBI and show 
how to generalize the path integral formalism to a many-body problem and subsequently compute the joint distribution 
of heights at equilibrium. In particular, for a large system, the joint distribution is shown to have the Wishart form 
in Eq. ([T]) with parameters /3 = 2 and TW — TV = a— 1/2 where a appears in the potential in Eq. ([2]). 

A. One interface 

Let us first consider the case of one single interface (TV = 1). The interface is described by its height h{x) for x 
from to L. When we think of the interface as a walker (or Brownian motion), the height h plays the role of the 
position of the walker, while the coordinate x along the substrate corresponds to time. The substrate can then be 
seen as a wall at height zero: one has h{x) > for every x. In the stationary state at thermal equilibrium, the energy 
of a configuration {h{x)} of the interface can be expressed as 

E[{h{x)}] = E.iast [{h{x)}] + U[{h{x)}] (3) 
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with Ef,iast [{^(2;)}] = \ Jq dx being the elastic energy (or the kinetic energy of the walker) and — 

dx V{h{x)) describes the potential energy due to the interaction potential V{h) with the substrate. The statistical 
weight of a configuration [{ft,(a;)}, < a; < L] of the interface is thus simply (setting ksT = 1 where ks is the 
Boltzmann constant and T the temperature) given by the Boltzmann weight 

P [{h{x)}] oc exp { - E[{h{x)}] } oc exp |-i (S''^'' ^ L^^ ^'^^''^^ "^^^ ^'^^ 

We assume periodic boundary conditions: h{0) = h{L) = h > 0. 

In absence of an external potential {V = 0), the interface is depinned, with a roughness exponent x = 1/2, which 
means (ft.^) — (ft.)^ oc L^^ cx L. In this case the interface is just the trajectory of a free one dimensional Brownian 
motion: the displacement h of the walker grows as the square root of the 'time' (i). But when the substrate induces 
an attractive potential, the interface remains pinned to the wall (substrate). The interface then becomes smooth 
(x — 0) in this case. 

Given the overall statistical weight of the full configuration of an interface in Eq. ^ over x € [0, L] , our task next is 
to compute the 'marginal' height distribution P{h) of the interface at a fixed point x in space, by integrating out the 
heights at other points. Note that due to the translational symmetry imposed by the periodic boundary condition, 
this marginal height distribution P{h) is independent of the point x, which we can conveniently choose to be x = for 
example. This integration of all other heights except at (or L) can be very conveniently carried out by the following 
path integral that allows us to write the marginal pdf P{h) as 

P{h)^ r^^^ "m(x)e-^[f'(^»] l,(,)>o (5) 

Jh(0)=h 

where the symbol lh(a;)>o is an indicator function that enforces the condition that the height at all points x G [0, L] 
is positive and the energy i?[{/i(a:)}] is given in Eq. (|4]). 

The path integral can be reinterpreted as a quantum propagator: 

P{h)o^{h\e-^^\h) (6) 

with the Hamiltonian 



iI=4;^ + VW (7) 



and with the constraint h > Q. The problem is now the one of a quantum particle in one dimension with position 
h{x) at time x, described by the Hamiltonian H (in imaginary time). 

We assume now that the energy spectrum of H is discrete (this will be the case in presence of the confining 
potential). The propagator can be decomposed in the eigenbasis of H: 

Z^E ^ 

where ipE is the eigenfunction of energy E. Thus calculating the marginal height distribution is equivalent, thanks 
to the relation in Eq. ([S]), to calculating the full spectral properties (i.e., all eigenvalues and eigenfunctions) of the 
quantum Hamiltonian H . In Eq. ([6]) and ([8]), the size L of the substrate (in the classical system of interfaces) plays 
the role of the inverse temperature in the associated quantum problem (but it has nothing to do with the temperature 
of the interfaces). Hence, in a large system L ^ 00, only the ground state (V'So = V'o with energy Eq) contributes 
to the sum in Eq. ([8j). Henceforth, we will always work in this limit where the marginal pdf is given by the exact 
formula 

P{h)^\Mh)\^ (9) 

Let us make a quick remark here. While the results in Eqs. ([5]) and Q may apriori look evident, they are however 
a bit more subtle. For example, the r.h.s. of Eq. ([9]) is, in a quantum mechanical sense, the probability density of 
finding a particle at h in the ground state. But it is not obvious (and needs to be proved as done above) that it also 
represents the height distribution of a classical model. 
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Thus our task is now to determine the exact ground state of the quantum Hamiltonian H. Analytically this is only 
possible for integrable H. With the choice of potential V{h) as in Eq. the quantum Hamiltonian H is fortunately 
integrable. The eigenfunction ^pn{h) satisfies the Schrodinger equation: 



1 rl"^ I 

H^n = -T^-^+V{h)ijn=En^n (10) 



with the boundary conditions, ?/'„(ft. = 0) = (due to the hard wall at /i = 0) and il)n{h oo) = 0. The solution is 
of the form 

Mh) =cne'^''' h" d"~^\bh''); with £;„ = 6(27i + a + i) (11) 

with n a non- negative integer (discrete spectrum), c„ a normalization constant and £« ^ a generalized Laguerre 
polynomial of degree n 

«-)=t(r-i)^ m 

i=0 ^ ' 

Note that for 7 = 0, (a;) = Ln{x) reduces to the ordinary Laguerre polynomial 

n \ (— a;) 



^»(-) = E UPtt^- (13) 



We also note that the generalized Laguerre polynomial in Eq. (|T2)) can alternately be expressed as a hypergeometric 
function 

>CI(x) = ^^^±^iFi(-n;7+l;x). (14) 
where (a)„ = (a)(a + l)-.-(a + n — 1) is the Pochhammer symbol and 

Ma,, a,; 6„ 5,; .) = g (5,)„...(5,)„ ^^^^ 
Finally, as 'Co(^) = 1' P'^'^ '^^ height of the interface is (in the limit of a large substrate L 00) 

P{h) ^ \Mh)\^ = \co\' e-"^' h'" (16) 

with|co|2- 



r(Q+i/2) ■ 

The mean m = (h) and variance = Vai{h) — {h^) — (h)^ of the interface height are thus easy to compute: 



= (h) = r dhh P{h) = J^'^^^^ (17) 

io ' ' V6 r(a + 1/2) 



2 l + 2a 1 / r(a + l) y , , 

ai = Var(/i) = — ^ (18) 

^ ^ ^ 2 6 6 Vr(a + l/2)y ^ ' 

B. interfaces 

Let us consider now N nonintersecting (1 + l)-dimensional interfaces over a substrate of size L. The z*'* interface 
is described by its height hi{x) for x from to L. Since the interfaces are nonintersecting, we can assume that they 
are ordered: < hi{x) < h2{x) < ... < h^^x) for every x. The only interaction between the interfaces is their 
fermionic repulsion (they do not cross). However, we will see that this constraint drastically changes the statistics of 
the interfaces. In the stationary state (at thermal equilibrium), the energy i?[{ft,(j;)}] of a configuration {h{x)} of 
one of the interfaces is given by ([3]) (same form as we assumed for one single interface), with an elastic energy and a 



8 



potential V{h) again given by Therefore the statistical weight of a configuration {hi{x) ;1 < i < N, 0<a;<L} 
of the whole system is simply (setting kgT — 1 for simplicity): 



P [{hi{x)}i^^] cx exp 



^-2 



(19) 



We assume again periodic boundary conditions: for every i, hi{Q) = hi{L) = hi with < /ii < /12 < ••• < /itv- The 
configuration space can thus be seen as a cylinder of radius L/2'k. 

The joint probability distribution of the heights of the interfaces at a given position (position x that can be taken 
to be by cylindrical symmetry, as we already noticed) can again be expressed as a path integral: 



P{hi,h2,...,hN)^T\ [ ^ 'm,(a:)e-^'^[f'-("»] l„„(,)>...>„,(,)>o (20) 



f.hi(L)=hi 
lhi{0) = hi 

The path integral can then be reinterpreted as a quantum propagator for N particles: 

P(/ll,/l2, ...,/lAr) OC (/ll,/l2, ...,/lAr|e"^^|/li,/l2, ...,/lAr) (21) 

where the many-body Hamiltonian is given by 

with ^ = E^^ = -E^J^+E^('^^) (22) 

i i i 

with the constraint > ... > hi > 0. The problem is now the one of N independent fermionic particles in one 
dimension with positions hi{x) at time x, described by the single particle Hamiltonian Hi (in imaginary time). 

Exactly as for one single interface, the propagator can be decomposed in the eigenbasis of H (Hamiltonian for N 
particles) and the joint distribution of heights is given by 

PihiM, ...,/..) = \Mh..h, h,)f 

where ?/'-e(^1: ^2, ■ • ■ , ^at) is the many-body wavefunction at energy E. Analogus to the single interface case, when 
the size of the system L tends to infinity, only the ground state (iV-body wavefunction) contributes to the sum. 
In this limit, the joint probability is simply: 

P{hiM.-,hN) = \^o{hi,...,hN)\'' (24) 

As in the single particle case, we emphasise that the relations in Eqs. ((23|) and ([24|) may apriori look evident, but 
they need to be proved as the l.h.s and r.h.s. of these equations refer to the probability density in a classical and 
a quantum problem respectively. We note that in the context of step edges on vicinal surfaces (in the absence of a 
wall), the relation p4|) was implicitly assumed in Ref. [3], but not proved. 

In the case of one interface, we computed the single particle wavefunction ipn in (llip . As the particles are indepen- 
dent fermions, the (iV-body) ground state wavefunction "^^is & N y. N Slater determinant. It is constructed from the 
N single particle wavcfunctions of lowest energy, the i/'i for i from to — 1 : 

*o(/ii,...,W)(xdet(^,_i(/ij))cxe-^2:.'^^. del (ct'^^'\bh])) (25) 

k 

Note that C^^~^\hh'^) is a polynomial of h^ of degree i — Any determinant involving polynomials can be reduced, 
via the linear combination of rows, to a Vandermonde determinant which can then be simply evaluated. We then get 

P{hi,...,hN)^\^o{hi,...,hN)? o^e-'^^'^i Ylhl" Ylihj-h',y (26) 

k i<j 

where /i^'s are positive. Note that due to the symmetry of the above expression, the ordering constraint hi < ... < hjq 
can be removed by simply dividing the normalization constant by A^!. 

For interfaces that are allowed to cross, the joint probability distribution of the heights has a similar form, but 
without the Vandermonde determinant: P(/ii, /iat) = P{hi)...P{his[) cx e^^'^^'^k J|j,/i^". The Vandermonde 
determinant Y\i^j(h1 — h'j)'^ comes from the fermionic repulsion between the interfaces. In particular, one has (as 
expected) for nonintersecting interfaces P{hi, ...^h^) = if /li = hj for i 7^ j. The consequence of this repulsion on 
the typical magnitude of the heights of interfaces will be explored in the next section. 
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C. Relation to the Eigenvalues of a Wishart Matrix 

We recall from the introduction that an {N x N) Wishart matrix is a product covariance matrix of the form W 

X where X is a Gaussian {M x N) rectangular matrix drawn from the distribution, P{X) cx exp —^Tt{X^X) 

For M > N, all eigenvalues of W are non-negative and are distributed via the joint pdf in Eq. ([T]). In the "Anti- 
Wishart" case, that is when M < N, W has M positive eigenvalues ( and N — M eigenvalues that are exactly zero ) 
and their joint probability distribution is simply obtained by exchanging M and N in the formula ([1]). 

The joint probability distribution of the heights of the interfaces in our model in Eq. can then be related to 
the Wishart pdf in Eq. ([1]) with f3 = 2 after a change of variables hhf = \i 

P{hi,...,hN)dhi...dhN (xe-''^"'''^ Y[hl" Y[{h^ - h^j)^dhi...dhN 

k i<j 

= Ajve-^-^'= J|a;^"^ l[{K-XjfdXi...dXN (27) 

k i<j 

where is a normalization constant. Recall that the parameter a > 1. By choosing a = (M — N) + i > 1, one 
recovers the Wishart pdf in Eq. ([T]) for (3 — 2 (complex matrices X) and arbitrary M — N > 1/2. Thus by tuning the 
amplitude a of the repulsive part of the potential in Eq. ^ one can generate the Wishart ensemble with a tunable 
M (with M — N > 1/2). The normalization constant An can be computed using Selberg's integrals [s^ and one gets 

A],' = n (fc! r(« - ^ + fc)) (28) 



Wishart pdf with arbitrary /3 > and M > N: We note that our model above generates a Wishart pdf with 
arbitrary M ~ N > 1/2 but with fixed (3 = 2. It is possible to generate the Wishart pdf with arbitrary (3>Q also by 
introducing an additional pairwise repulsive potential between the interfaces. For example, we may add to the energy 
functional in Eq. (jl9p an additional pairwise interaction term of the form, 



E 



l<j<k<N 



Vpair {hj{x), hk{x)) dx 



(29) 



where the pair potential V^air(^j, h^) has a specific form 



Vpair(^j, /ife) 



(3 f(3 



(30) 



with (3 >Q. In this case, once again using the path integral formalism developed above, we can map the computation 
of the joint distribution of heights to calculating the spectral properties of a quantum Hamiltonian via Eqs. (j23|) 
and p4|) . The corresponding quantum Hamiltonian turns out to be exactly the Calogcro-Moser model [36] which 
is integrable [stI . |38| . In particular, using the exact ground state wavefunction of this Hamiltonian we get our 
corresponding joint distribution of interface heights in the following form 



p(;ii,...,;ijv) = |*o(/»i,...,/»Ar)p(xe-''^'=''' \{hl^ \{{h 



(31) 



which, after the usual change of variables bh^ = Xi, corresponds to the general Wishart pdf in Eq. ([T]) with arbitrary 
/3 > and a tunable M — N = {2a + 1 — (3)/(3. We note that the procedure used above to obtain a variable (3 
random matrix ensemble was used before in the context of step edges on vicinal surfaces without a hard wall where a 
corresponding Gaussian matrix ensemble with tunable (3 was obtained Q . 

Note that for interfaces in presence of a wall, while the first term in the pair potential in Eq. pop is quite natural 
and can arise out of entropic origin as well as dipolar interaction between step edges Q , the second term however does 
not have any physical origin. Unfortunately if one gets rid of this term, the integrability of the quantum Hamiltonian 
also gets lost. In any case, in the following we would focus only on the physical /? = 2 case. 
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III. STATISTICAL PROPERTIES OF THE MODEL 



Once the joint distribution of heights at equiHbrium is known, one can, at least in principle, compute the statistics 
of various relevant quantitites such as the average density of lines at height h, the distribution of the maximal and the 
minimal height, the distribution of the center of mass of the interfaces etc. In this section we show how to carry out 
this procedure and derive some explicit results upon borrowing the techniques developed in the context of random 
matrix theory. Many of the results in Sees. Ill A-C actually follow from a simple change of variables in the already 
known results of random matrix theory. For the sake of completeness, we remind the readers some of these results 
from random matrix theory and draw the consequences for our interface model. However, Sec. Ill D (where we derive 
the distribution of the center of mass of the interfaces) presents completely new results. 

For simplicity, we will focus here only on large system (L —^ oo) properties. In particular, our focus would be to 
understand the effect of fermionic repulsion between the interfaces (nonintersecting constraint) and also the effect of 
the external confining potential on the statistics of the above mentioned physically relevant quantities. 

We have shown in the previous section that the joint pdf of interface heights (hi), after the change of variables 
bh^ = Xi, is the same as the joint pdf of Wishart eigenvalues (A^) in Eq. ([1]), which can be re- written as a Boltzmann 
weight 

Pn{Xi, Xn) oc exp [-l3E^s{Xi, Ajv)] (32) 

with the effective energy 

^^off - ^ ^ Afc - a ^ In Afc - ^ In | A, - A, | (33) 

k k i<j 

where a = ^l+MiiiY. — In this form, the A^'s can be interpreted as the positions of charges repelling each other 

via the 2-d Coulomb interaction (logarithmic), but are confined on the 1-d positive axis and in presence of an external 
linear -|-logarithmic potential. The Dyson index P plays the role of inverse temperature. Our model of interfaces 
corresponds to /3 = 2 (see Eq. ([27]) V 

We already noticed that the nonintersection constraint for the interfaces is equivalent to the presence of the Vander- 
monde determinant, and thus the logarithmic Coulomb repulsion, in the joint probability distribution. For independent 
interfaces (allowed to cross), there is no Vandermonde term and hence the logarithmic repulsion term in Eq. ([55)1 is 
absent. In that case, balancing the first two terms of the energy gives a typical height of order one: /ityp ^ 0{1). But 
for the nonintersecting case, when the number TV of interfaces becomes large, the logarithmic repulsion is stronger 
than the logarithmic part of the external potential (provided a is not proportional to N). Therefore, balancing the first 
and the third term in the effective energy gives, for large N, iVAtyp ^ N'^, thus Atyp ~ iV or equivalently /ityp ^ ^/N. 
The effect of repulsion is strong: the interfaces spread out considerably. 

Below we first compute the average density of states, followed by the computation of the distribution of the topmost 
interface (maximal height) and the lowest interface (minimal height) that is the closest to the substrate. Finally we 
analyze the distribution of the center of mass of the heights. 



A. Average density of states 

We would first like to know what fraction of N interfaces lie, on an average, within a small interval of heights 
[h, h + dh]. This is given by the average density of states (normalized to unity) 

1 ^ 

PNih)^^Y.(^{h-h,)) (34) 

4=1 

As we explained above, we expect the typical height scale to be of order htyp ^ Vn for a large number N of interfaces. 
Furthermore, the density of states is normalized to unity: dh pN{h) — 1. Therefore the density is expected to 
have the following scaling form for large N: 

Our goal is to compute this scaling function g{x). This can actually be simply read off from the known results on 
Wishart matrices which we now recall. Consider the Wishart matrix with M > N with eigenvalues distributed via 
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the joint pdf in Eq. ([T]). In the asymptotic limit N oo, M ^ oo keeping the ratio c = N/M fixed (with c < 1), the 
average density of states of the eigenvalues is known [s^ to be of the form: 

where the Marcenko-Pastur scaling function f{x) depends on c (but is independent of /?) 



fix) = —^(x+-x)ix^x-) (37) 

which has a non-zero support over the interval x G [a;_,x+] where x± — (^-^ ± 1^ . Note that in the limit c 1, 

which happens when M ~ N ^ 0(1) for large iV, x_ — > and x+ 4. 

Our interface model, after the customary change of variable bh^ — Xi, corresponds to the Wishart ensemble in Eq. 
^ with /3 — 2 and M — N = a — 1/2. Hence, as long as a ~ 0(1) for large N, c = N/M — 1 in our model. Using 
bh"^ = A and c = 1 in Eq. ^ and the average density of states in the interface model then indeed has the 

scaling form in Eq. psp for large N with the scaling function 



,(.)^V^-.^ (38) 



where b is the frequency of the harmonic part of the potential (see The average density is a quarter of ellipse, 

as shown in figure [H It has a finite support , 2 ■ 

Thus the interface heights spread out for large A'^ as a result of the nonintersection constraint. Let us compare this 
result to the case of independent interfaces that are allowed to cross each other. In that case, the average density 
of states is simply pN{h) = P{h) cx e~^^ h?°' . It is independent of N (evidently!) and has a non-zero support 
over the whole positive h axis. It vanishes when h tends to zero, and rapidly decreases to zero when h becomes 
large. Thus most of the interfaces lie on an average close to the wall at a distance of 0(1). In contrast, the heights 
of nonintersecting interfaces, on an average for large N , have a compact support over a wide region. The density 
vanishes at the upper edge as a square root singularity and the upper edge itself grows as -s/iV, thus spreading the 
interfaces further and further away from the wall as N increases. 

The average of the upper height (maximum) is expected to be given by the upper bound of the density support: 

{hmax) ~ 2 \J^- The lower bound of the support of the density is zero for large N in first approximation. We will 
show more precisely that the average height of the lower interface (minimum) is proportional to 



bN 

Finally, the average of all interface heights (h) — ( (^'i+''2+-..+^tiv) ^ ^^^^ computed for large N- 







(h)^ I hpN{h)dh^^^j (39) 



This differs drastically from the case of independent interfaces, where the average of all heights is the same as that 
for one single interface: (h) = ^^-^'{^+1/2) ^^^^ section III Ap is independent of N, but depends on a (parameter 
associated to the part of the potential proportional to ^). In contrast, for nonintersecting interfaces, the part of 
the potential proportional to -j^ becomes negligible compared to the repulsion between interfaces, thus {h) does not 
depend on a, but grows with N (see (|39)) ). 



B. Maximal height of the interfaces 



In this subsection we compute the distribution of the height of the topmost interface (maximal height), the one 
furthest from the substrate. The average of the maximal height is given by the upper bound of the density support 
(see section UlI A^ : 

(/imax) « 2 7^ for large N. (40) 
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But we would like to know the full distribution of the height of the topmost interface, not just its average. For that 
purpose, we can again take advantage of the mapping between our interface model and the Wishart random matrix. 
Under this mapping, the height of the topmost interface /imax is related, via the change of variable bh^^^ = Amax, to 
the largest eigenvalue A^ax of the Wishart matrix. The distribution of A,„ax has been studied in great detail and we 
can then directly use these results for our purpose. 

Let us recall briefly the known properties of the largest eigenvalue Amax of Wishart matrices whose eigenvalues are 
distributed via the pdf in Eq. ([T]). For our purpose we will only focus on f3 ~ 2 and M > N with M — N ~ 0(1). 
In this case the parameter c = N/M tends to the limiting value c = 1 for large N, indicating that the upper edge 
of the Marcenko-Pastur sea, describing the average density of states in Eq. ([37|) , approaches x+ ^ 4 and lower edge 
X- 0. Thus, the average of the maximal eigenvalue of a Wishart matrix is (Amax) ~ 4A^. Furthermore, Johansson 
[s^ l and Johnstone independently showed that the typical fluctuations of Amax around its mean iN are of order 
i.e., 



A„ 



(41) 



IS 



where the random variable X2 has an iV-independent distribution for large N, Prob(x2 ^ x) = ^2(2;) where F2{x 
the celebrated Tracy- Widom distribution |3J] for (3 = 2. However, for finite but large N , the tails of the pdf of Amax 
(for |Amax — 4A^| ^ 0{N)) show significant deviations from the Tracy- Widom behavior. The behavior in the tails of 
the pdf P(Amax — t, N) is instead well described by the following functional forms [3^ . valid for arbitrary /3, 



P{t,N) ~ exp 



exp 



2^ f^N-t 



V N 

t-AN 
N 



for t < 4iV ; 
for i > 4iV ; 



(42) 
(43) 



where ^±{x) are the right (left) large deviation (rate) functions for the large positive (negative) fluctuations of Amax- 
Interestingly, the explicit form of the rate functions $± {x) have recently been computed and were shown to be 
indepedent of /3. The left rate function <I> jx) was computed in Ref. [40| using a Coulomb gas method developed in 
the context of Gaussian random matrices |4l| and is given for a; > by 



$_(x) = In 



X X 

8 ~ 64' 



The right rate function <I>_|_ (x) was also computed very recently using a different method 



<f+(x) 



hi{x + 4) + 



1 



G 



(44) 



(45) 



where G{z) =3 F2 [{1, 1, 3/2}, {2, 3}, z] is a hypergeometric function. For small argument x, the two rate functions 
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have the following behavior [40|, |42 



^^(x) w x3/384 



(46) 
(47) 



Using these results, it was shown [40, |42| that both large deviation tails of the pdf of Amax in Eqs. (|42|) and (|43|) 
match smoothly with the inner Tracy- Widom form. 

These results can then be directly translated to our problem of interfaces identifying bhf-^^^^ — Amax- For large N, 
the typical fluctuations of /imax around its mean are Tracy- Widom distributed (see figureS]) over a scale ^ 0{N^^^^). 
More precisely, we get 



N 



2 3 b 2 TV 6 



(48) 



where X2 is Tracy- Widom distributed, Prob(x2 l£ x) — F2{x). 

This gives in particular the first finite size correction to the leading term for the average of the maximal height, in 
the large N limit: 



2a/ y +2-3 (X2) 



(49) 



where {X2) ~ —1.7711 [3J|. The variance can also be computed from 
distribution (xi) - (X2)^ « 0.8132 [H: 



and the known variance of the Tracy- Widom 



Var(/i 

n 



(^max) 

2-4/3 
b 

0.32 

5^1/3 



{K 



(50) 



PN{h) 



Tracy- Widom 




FIG. 4: The dashed line shows the density of states pN{h) (quarter of ellipse). The height of the upper interface (maximal 
height) is centered around its mean 2-J ^ with fluctuations of order 0{N^^^^) described by the Tracy- Widom law. 



Similarly, the atypical large fluctuations of /imax around its mean (for |/imax — '^\^N/b\ ^ 0(iV^/^)) are described 
by the large deviation tails as for A^ax in Eqs. (|^^ and with /? = 2. For the left large deviation, we get for 

6^2 _ 4iv ~ 0{N) with bt^ < m : 

P[Vax^^,iV]^exp|-2jV^$_ ( ^^^^^' )| (51) 



where $_(x) is given in Eq. 
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P (/imax = 2xVN 
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FIG. 5: Pdf of the maximum of the heights (large deviation), plotted for b — 1 and for different A'^, as a function of the rescaled 
height X = As increases, the rescaled location of the peak of the pdf approaches to x ^ 1. The width of the regime 

around the peak where the Tracy- Widom law is valid is reduced to ~ jsf-''-/^-''-/^ ^j^jg gg^le. The rest of the pdf 

beyond the peak are described by the large deviation tails in Eqs. (|5ip and (|52p . 



Analogously the large and rare fluctuations to the right of the mean can also be computed from the exact expression 
of <&+(a;) in Eq. ((i5l) . Replacing again A^ax by bh"^^^, we get the right large deviation tail of the pdf of /imax, for 
large N and for bt^ - 4iV - 0{N) with bt^ > 47V: 

P[/jmax = i,Af] ~cxp|-27V$+ f (52) 



where $+(a;) is given in Eq. (|¥5l 



C. Minimal height of the interfaces 

We have seen in subsection IIII Al that the lower bound of the support of the average density of interfaces is indeed 
zero in the first approximation as iV — > oo. Since the lower edge of the support is also precisely the average height of 
the lowest (close to the substrate) interface, we have (ft-min) — > as — > cxd. This is clearly an effect of the fermionic 
repulsion between the interfaces, because for 'independent' interfaces (that are allowed to cross) the height of the 
lowest interface (minimal height) is of order /imin ^ 0{1): it does not see the other interfaces. The hard wall at the 
origin in the problem of interfaces corresponds to the constraint that the Wishart eigenvalues must be positive. In 
the Wishart ensemble, the neighborhood of the origin is, thus, called the hard edge of the spectrum. The distribution 
for the hard edge is known to be related to the Bessel kernel [H, |43|, just like the distribution at the soft edge (for 
the maximum eigenvalue) is related to the Airy kernel. 

To know more precisely how Zimin decreases with increasing A^ in presence of the 'nonintersection' constraint, we 
need to find the statistics of /imin for large but finite N, which is precisely the objective of this subsection. The 
main result of this subsection is to show that for nonintersecting interfaces, the minimal height is typically of order 

hmin ~ O to leading order in large A^. For special values of the parameters, we are also able to calculate the 

full distribution of the minimal height /imin as discussed below. 

Under the customary change of variables = , it follows that /imin has the same distribution as \/ Amin/^ where 
Amin is the minimum eigenvalue of Wishart ensemble with parameters (3 = 2 and arbitrary M > N. The minimum 
eigenvalue of the Wishart ensemble has been studied before p[5| . with applications in the quantum entanglement 
problem in bipartite systems [i^. When M — N ^ 0(1), which is precisely our case since a = M — N + 1/2 ~ 0(1), 
the minimum eigenvalue is known to scale, for large A^, as Amin ~ 1/A^ for arbitrary /3, though it has only been proved 
exactly for special values oi M — N and j3, e.g., for /3 = 1 and M = A^ or for /3 = 2 and M = N [i^. For the interface 
model it then follows quite generally that /imin ~ I/a/a" for large A^ for general a ~ 0(1). The fermionic repulsion 
between interfaces has thus again a strong effect: the lowest interface is pushed very close to the substrate since its 
height is of order h^i^ ~ I/n/a" for large A^ (instead of 0(1) for non-interacting interfaces that can cross). 
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To go beyond this scaling behavior for large N and compute precisely the statistics of /imin for arbitrary A'' seems 
difficult for general a. Below we show that for the special case a — 3/2, it is possible to compute the full distribution 
of /imin for all N. 

It turns out to be convenient to compute the cumulative distribution function (cdf) of the minimal height 
Prob [/imin > VCj ^] J for arbitrary (^. We use the notation for the convenience of scaling as seen below. Our 
starting point is the central result for the joint pdf of interface heights in Eq. (pS)) . Clearly, the event that the 
minimum height /imin > VC is equivalent to the event that all the heights are greater than y^: hi > ^JQ for all 
i = 1, 2, . . . , A^. Hence, 



Prob 



dhi 



dhN P{hi, /i2, ■ • ■ , hN) 



(53) 



where the joint pdf P{hi,h2, 
have 



, hjsi) is given in Eq. (|26p . Making the standard change of variables, hKf — we then 



Prob 



A 



N 



(54) 



where the normalization constant is given in Eq. 
(|54p in a more compact form 



k i<j 

Next, making a shift A^ — bC, + Xi, one can rewrite Eq. 



Prob 



= Aat e 



-bNC . 



where the function w{z) is given by the multiple integral 



w(z) 



/•oo noo ^ 

/ dxi... dxN e-'^-'^" J|(a;fe + z)"-3 J|(x, - x 

•^0 -^0 fc=l ^<J 



(55) 



(56) 



For notational simplicity we have suppressed the N and a dependence of w{z). The multiple integral in Eq. (|56p is 
not easy to evaluate for general values of the parameter a. However, one can make progress for special values of a. 

We first note that when the parameter a— ^=M — TV is an integer, w{z) is a polynomial of z of degree a — ^. In 
the special case a = 3/2, i.e., M = + 1, one can explicitly evaluate w{z) by following a method similar to the one 
used by Edelman [4^ to compute the distribution of Amin for Wishart matrices with (3^1 and M = N. In this special 
case M — N + 1 and (3 = 2, the distribution of the minimum eigenvalue, in the large- limit, was already computed 
by Forrester [43l |. However, in this special case, one can actually calculate the distribution and the moments even for 
all finite A^ as we demonstrate below. For a — 3/2 and (3 = 2, we first compute two derivatives w'(z) and w"(z) of 
the function w{z) in Eq. (|56p . Using integration by parts and some rearrangements, we find that w{z) satisfies an 
ordinary second order differential equation for any A^ 

zw"{z) + {1 + z)w'{z) ~ Nw{z) =0 (57) 

whose unique (up to a constant) solution is in fact the ordinary Laguerre polynomial with negative argument 

N 



/ N \ z'' 
wiz)c^CN{-z)^Y.( k ] u 



(58) 



and finally, since Prob [/imin > 0, A'] = 1, we get 



Prob [/i„,i„ > t, iV] 



CNi-bt^) = e 



-bNt^ 



N 

E 

fe=0 



fjk ^2k 
k\ ■ 



(59) 



Hence, for the special case a = 3/2 we can then give an explicit expression for the pdf of /imin valid for all A^, 



P {hnun ^t,N) = -^Prob [/i,ni„ >t,N]=2b' f 



2 ,3 „~bNt^ /.(2) 



(60) 



where £^ (x) is the generalized Laguerre polynomial already defined in Eq. ([T 
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From the exact pdf in Eq. (l60l) one can calculate all its moments explicitly as well (see appendix- A for details). 
Wc find for the fc-th moment, for arbitrary N, 

(hU = ^^^^ 2F, {-{N - 1), k/2 + 2; 3; -1/A^) . (61) 

One can then work out the asymptotic behavior of the moments for large N. For example, one can show (see 
appendix- A) that the average value (fc = 1) (/imin) decays for large N as 

« (62) 

with the constant prefactor ci given exactly by 

/ TTG 

ci = W — /o(l/2) = 1.5538... (63) 
where /„(z) is the modified Bessel function of the first kind with index n 

1 / 2'\ 2k+n 



fc=0 



For large N, one can also work out precisely the scaling behavior of the full pdf of /imin given in Eq. ([50]) and 
recover the result of Forrester Since typically hmin ^ 1/v^i one expects that its pdf (normalized to unity) has 
a scaling form for large N 



P (Vin = t,N)^ VNUi,,{t VN). (65) 
The scaling function /min(a;) can be computed explicitly from Eq. (|60p . We get 

/„;„(.) = ^lun^ -^P U^,^ ^7^^^ 



= lim b'^x^e-"^" ,Fi(l-7V;3;^) 

= 2b xe-'"'' h{2xVb). (66) 



This function has the following asymptotic behavior 

fmin{x) K b^ x^ as X ^ 



63/4 

A plot of this scaling function is given in Fig. 6. 



f^^^bx^+2^x asx^oo (67) 



D. Center of mass 



We study in this subsection the distribution of the center of mass of the heights 

^ h±_±]^ (68) 

for large N . If the interfaces were allowed to cross, the heights hi would be independent and identically distributed 
variables (i.i.d.). In that case, the distribution of the center of mass Gat would be a pure Gaussian distribution in the 
large N limit (central limit theorem). But in our model, due to the repulsion between the interfaces, the interface 
heights are strongly correlated. What is the effect of the repulsion (nonintersecting constraint) on the center of mass? 
We will show how to compute the pdf of the center of mass for large N upon borrowing some techniques developed in 
the context of random matrix theory, using in particular the analogy between the Wishart eigenvalues and a Coulomb 
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FIG. 6: Rescaled probability density function of the minimum of the heights, in the case a = 3/2, plotted for & = 1 and for 
N = 5 (dashed line) and A'' = 20 (dotted line) and compared with the limiting distribution /min(a;) (solid line). In the large A'' 
limit, indeed the curves approach the limiting distribution fmini^)- 



gas of charges. We will show that the pdf of the center of mass P{Gn = i^) has an extraordinarily weak non-analytic 
behavior at i/ = {Gn) (where (Gn) is the average of the center of mass), which is shown to be a direct consequence 
of a phase transition in the associated Coulomb gas problem. 

Since the typical height of an interface ft-typ ~ for large N, it follows that the center of mass Gm « 0{y/N). 
More precisely, by symmetry of the joint pdf of the heights, the average of the center of mass is given by the average 
height (see (p^ ): 

8 /~/V~ 8 

{Gn) = (h) ^ —\ —= fiVN for large N, where ^ = (69) 

ovr V 3nvb 



Let us thus write v — svN, where the scaled variable s ^ 0(1). 

The main result of this subsection is to show that in the scaling limit iV — s- oo, — > cx) but keeping the ratio 
s = j^/Vn fixed, the pdf of the center of mass scales as: 



P {Gn — v) (x. exp 



V 

~N 



(70) 



where the associated large deviation function $(s) is plotted in figure [U has the following asymptotic behavior 



-21ns for s ^ 0+ 

$(s) « <; (71) 

h for s +CX) 

and is a non-analytic smooth function: <l>(s) is infinitely differentiable everywhere but it is not analytic. More precisely, 
we will show that $(s) is given by 

(i>(s) = I ^+1^1 s > where $~ and <i>^ are analytic functions on their domain of definition^ (72) 

andwith $+(s) - $-(s) « -7rV6(s-/^)e"-^^(-f) e^^'^^-i) as s ^ /x+ (73) 
$(s) has thus an essential singularity at s = it is not analytic. 



$ is analytic on the complex plane except for a branch cut along the negative real axis; and <!>+ is analytic on an open set of the 
complex plane including the (real) half-line s > fj,. 
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But all the derivatives of $ exist and are continuous. In particular, for s — > /i, $ has, in first approximation, a 
quadratic behavior: 



(s - liY 

'^{s) w ^ ^ for s — > /X with /.t 



— and cF=-\ \ 



(74) 



The pdf of the center of mass can thus be approximated by a Gaussian around its minimum (s = /i), which gives the 
mean and variance of the center of mass: 



{G.)^,^^l^l^ and VVa^ ^ ^ {Gl) - {G.? ^^^\\ff, 



3n\l b 

We will also derive an exact closed form for (f>^(s): 



(75) 



^ ' 32 



L(s] 



1 



where L{s) 



2^^^ 5i(s)"^^^ -21/^51(5)1/3 with gi{s) = -Sirs Vb+y/l6 + 9 b 



But we will see that $+(s) is more difficult to compute: we will only derive its asymptotics (s 



and 
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To derive these results, let us start with the pdf of the center of mass: 
P{Gn^sVn)^ f dhi... f dhN s( '^^^ "^'^^''^ - sVn\ P{hi,...,hN) 



^0 



= rdx,... rdxj^e-^^'^ nxr- ^(^..-A,-)^^f ^+^•^+^ -.^/jv) (77) 

JQ Jo ^ ~g \ Nvb J 

where we have used Eq. ((27)) . The integrand (without the delta function) can be written as exp [— £^{Ai}] where 

N / 1\ ^ 

i?{AJ = ^A.- U--^ln(AO-^ln|A, -A,| (78) 

i=l ^ ' i=l j^k 

can be interpreted as the energy of a Coulomb gas of N charges with coordinates {Xi} as mentioned earlier in section 
III. Thus the calculation of the distribution of the center of mass reduces to the calculation of the distribution of 
a particular functional of this Coulomb gas. This can be performed exactly for large N using a functional integral 
method followed by saddle point calculations. This method has been used recently in several contexts: for example, to 
calculate the large fluctuations of the maximum eigenvalue of both Gaussian and Wishart random matrices 40, 4ll.[43|. 
to compute the purity distribution in bipartite entanglement of a random pure state [47l| and also to compute the 
distributions of conductance and shot noise for ballistic transport in a chaotic cavity |48j . 

To evaluate the multiple integral in Eq. (|77p by the functional integral method one proceeds in two steps. First step 
is a coarse-graining procedure that sums over (partial tracing) all microscopic configurations of {A^} compatible with 

a fixed normalized (to unity) charge density pn{X) = J2i H'*^ ~ ^i) and a fixed value of s = J2i V^/ 

The next step is to integrate over all possible normalized charge densities with fixed s-this is the functional integration 

which is then carried out using the method of steepest descent for large N. 

To proceed, we first scale the positions of charges, x = ^ such that x ~ 0(1) and define the charge density in 
the X space p{x) = S [x — j^). With this scaling, it is easy to check that while the first and the third term in 

the energy expression in Eq. ([75]) are both of order ~ 0{N^), the second term multiplying {a — 1/2) (corresponding 
to the external logarithmic potential) is of order ~ 0{N), as long as {a — 1/2) is of order ^ 0(1). Thus this term 
becomes negligible for large N for any A^-independent a and to leading order in large N, the a-dependence just drops 
out. Then, the coarse-graining procedure gives, to leading order for large N, 

P{Gn = sVn) cx Jv[p] e-^' EAp]+oiN) (79^ 
where the effective energy functional is given by: 

poo poo poo / poo \ / poo \ 

Es[p]^ J xp{x)dx-J J p{x) p{x') In \x - x'\dxdx' + r( J ^/xp{x)dx - s Vb) + D (j p{x)dx ~ Ij (80) 
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We have introduced two Lagrange multipliers, R and D, in order to take into account two constraints. The first 
(associated to R) enforces the condition Gn = ^^^'y^^ ^"^ = s\/N or equivalently ^/xp{x)dx = s^/h (it replaces 

the delta function in the expression of P {Gn ~ s\/n'^). The second (associated to D) enforces the normalization of 

the density p: J^^ p{x)dx — 1. 

The functional integral ([79|) is carried out in the large N limit by the method of steepest descent. Hence: 

P{Gn = sVn) oc exp [-N^ [pc]] (81) 

where pc{x) minimizes the effective energy: ^^^^'^^^^^ ~ 0. The saddle point density Pc(^) is thus given by the equation: 

x + RVx + D = 2 pc{x')ln\x-x'\dx' (82) 
Jo 

Differentiating once with respect to x leads to the integral equation: 

1 + JL^2V r P^dx' = 2 i/, [pj (83) 
zy X Jrt X — X 







Hx [Pc] is called the semi- infinite Hilbert transform of pc (and P denotes the principal value). It is not easy to invert 
it directly. However, the finite Hilbert transform [y] = 'P can be inverted using a theorem proved by 

Tricomi ^49j]. According to Tricomi, the solution of the integral equation 

f(x) = V [ ^^^dt with a<x<b, lal + 161 < oo (84) 

Ja 

is given by 



yix) 



'n^\Jx — a\/h — X 



t — x 



for a <x <b (85) 



where Cq is an arbitrary constant. Tricomi showed that Co then satisfies: tt ^^y{t)dt = Co. Wc will hereafter assume 
that the saddle point density pc has a finite support and use Tricomi's result. 

So, the steps we need to carry out are (i) to find the solution Pc{x) of the integral equation (I83p which will 
contain yet unknown Lagrange multipliers R and D (ii) fix R and D from the two conditions: Pc{x) dx — 1 and 
/o°° \/x Pc{x) dx — s^/b for a fixed given s and (iii) evaluate the saddle point energy Es[p^ which is then precisely (up 
to an additive constant) the large deviation function <i>(s) announced in Eq. ((70|) . 

Physically, as the effective (external) potential for the charges is of the form Vf{x) ^ x + R^ + D (see equations 
((80)) and ((82)) ). we expect a different behavior of the charge density Pc{x) depending on the sign of the Lagrange 
multiplier R. 

• For i? > 0, the effective potential Vf{x) = x + R^/x + Z) is an increasing function of x for x > Q with minimum 
at a; = 0. In this case, the charges will be confined near the origin. Therefore the density must be large for small x, 
decreasing as x increases and finally vanishing at a certain x = L. We thus assume that Pc{x) has a finite support 
over ]0, L] where L is fixed by demanding that the density vanishes at a: = L: Pc{L) = 0. 

• However, for i? < 0, the effective potential is minimal for x = xq = ^ > Q. The density must be larger around 
X — xq. In that case, pdx) will have a finite support over [Li, L2] with Li > and where Li and L2 are fixed by the 
constraints pc{Li) = = pc{L2). 

We will see later that i? > corresponds to the left side of the mean of the center of mass (s < /i), and i? < 
corresponds to its right side (s > p). Thus there is a phase transition in this Coulomb gas problem as one tunes s 
through s = p OT equivalently R through the critical value R = 0. The optimal charge density has different behaviors 
for i? > and R < 0. When expressed as a function of s, this leads to non-analytic behavior of the saddle point 
energy, i.e., the large deviation function <&(s) at its minimum s = p. 
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1. Case R>0 (s < fi) 

Let us begin with the case R > 0, that wiU be shown to correspond to s < fi (left side of the mean of the center 
of mass). In this case, the effective potential is minimal for a; = 0. We can thus assume that Pc has a finite support 
over ]0,L] where L is fixed by the constraint Pc{L) = 0. In this subsection, we compute the saddle point density 
Pc{x) and derive an exact closed form for the energy Es[pc\ (and thus the function ^(s)). From this explicit form, 
we work out the asymptotic behavior of <&(s) for s —^ and for s p^ . For s — *■ we will see that the pdf can be 
approximated by a Gaussian -and this will give the mean and variance of the pdf of the center of mass. 

The (normalized) solution Pc{x), with support over ]0, L], of the integral equation (|83p can then be obtained using 
Tricomi's theorem in Eq. (|85p . The resulting integral can be performed using the Mathematica and we get 



p.(x) = i-y^^ + ^argth(yr|) forO<.<L (86) 

where argth is the inverse hyperbolic tangent. 

As the density pdx) must be positive for all x £]0,L[ (it is a density of states, of charges), such a solution (with 
support over ]0, i]) can exist only for R> 0. For i? 7^ 0, we have indeed pdx) ~ ^ as x ^ 0"*". Therefore R 
must be positive: R> 0. Conversely, it is not difficult to see that for R> 0, the density given in Eq. is positive 
for all X e]0, i[. In this phase {R> 0), as figure [7] shows, the Coulomb charges are confined close to the origin: the 
interfaces are bound to the substrate. 
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FIG. 7: Density of states pc(x) (density of charges) of the Coulomb gas associated to the computation of the pdf P{Gn = s\/iV) 
of the center of mass, in the case s < p = (^ > 0), plotted for different values of s (and for 6=1). The effective potential 

seen by the charges is minimal for x — 0, thus the density has a finite support over ]0, L] and diverges at the origin. 

• When s tends to p 0.85 for b — 1 (i.e. the center of mass tends to its mean value), L tends to 4 and pc tends to the average 
value of the density of states {R —> 0). 

• When s < p and s decreases (i.e. the center of mass is smaller than its mean and decreases) , L < 4 and L decreases also: the 
Coulomb gas of charges is more and more compressed, the charges are more and more confined close to the origin. 



We want to compute the pdf P{Gn — s V N). The basic variable is thus s. There are also three unknown parameters: 
R and D are two Lagrange multipliers and L is the upper bound of the density support. These parameters will be 
determined by enforcing the three constraints pc{x)dx — 1, Jp°° y/xpc{x)dx = s Vb and Pc{L) = 0. 

Hence, the parameters L and R are solutions of the two following equations: 

^ + - = .V5 and R-^^-^ m 
These equations can be solved exactly. In particular, we obtain the following expression for L = L(s): 

i(s) = (-.91 (5)1/3 21/3 + 2^/3 5i(s)-i/3^^ with gi{s) = ^3TTsVb+ ^/TeTob^ (88) 
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The saddle point energy can then be computed (from equation ([80l) and using ([82)1 for the calculation of the Lagrange 
multiplier D) as a function of L = L(s): 

Finally the distribution of the center of mass, in the large N limit, is simply given by the steepest descent method 
P{Gn = sVN) oc exp [-N^ [pj] : 



P{Gn — v) oi exp 



L{sY ^ (L{s)\ 1 



with <i>(s) = -l^-ln( ^ 1 (90) 



with L = L{s) given in Eq. (|88p . The additive constant has been chosen for convenience such that the minimum of 
$ is 0. $ is thus a positive function. $(s) is plotted in Fig. [H As expected, the minimum of <i>(s) is reached for 
s = ^ , where Gn = fJ- VN = (Gn) -the average of the center of mass. 

Validity of the regime where the density has a support over ]0, i]: i? > 0, s < /i 

As we noticed above, the density pc must be positive for every < a; < i, which is equivalent to demanding 
that R > 0. And from Eq. (|57|) . one can easily show that the constraint i? > is equivalent to s < /x. Thus the 
expression of $(s) given in Eq. ([50)1 is only valid on the left side of the mean of the center of mass: z/ < ^\/N (or s < fi). 

Limit s ^ /i~ (i? ^ 0+ ): Gaussian approximation of the pdf 

For s — > ^~ , $ can be expanded about its minimum: 

^(s) ~ i^l"^ where n = — and a ^ - \ r (91) 
In this limit, the pdf of the center of mass can be approximated by a Gaussian: 

P{Gn = sVN) oz e ^ as s ^ fi" (92) 

For large iV, only the vicinity of s = p,, where $ is minimum, will contribute. Therefore, the Gaussian approximation 
above gives the mean value of the center of mass and its variance: 

(G^) ^ (h) ^ ^ (93) 



VVa^ ^{G%) - (G^)2 « « i^T (94) 

This differs again strongly from the case of independent interfaces. For interfaces that are allowed to cross (they are 
thus completely independent), the average of the center of mass (Gat) ~ (h) = m is of order one, and its variance is 
given by \/Var(GAr) = where m (resp. ai) is the mean (resp. variance) of one single interface (see section ing . 

Both m and cti depend on a and b: they depend on the whole form of the potential V{h) = ^ — h ^^^7^3^- But for 
nonintersecting interfaces, only the harmonic part of the potential (with frequency b) has a non-negligible effect for 
large N (the a-dependence drops out, as we explained at the beginning of the section). And the relative standard 

deviation is of order O (^:^) for independent interfaces against O (-^) for nonintersecting interfaces. The 

relative fluctuations are strongly reduced by the fermionic repulsion. 

Limit s -^0+ {R ^ +00) 

For ,s 0+, the upper bound L(s) of the density support tends to zero like s^: L{s) w tt^ 6 + O(s^) as s ^ 0+ 
and thus $ tends to infinity : 

$(s) w -2 In s - i - In (^'] + 0{s In s) as s ^ 0+ (95) 
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The probability density function thus tends to zero as a power law: 

P{G = sVn) cx s^^' as s ^ 0^ 



(96) 



To summarize, for s < /x = 8/^37r-\/6^, the large deviation function <i>(s) = $ (s) characterizing the form of the 
pdf of the center of mass Gat to the left of its mean value is given by Eqs. (190]) and ([88)l . and is plotted in Fig. [H 




M s 

FIG. 8: Large deviation function $(s) of the pdf of the center of mass, such that P{Gn = sViV) oc e"^ The minimum 

of $(s) occurs at s = /I = 8/ (^Siry/bJ which corresponds to the average value of the center of mass. We have chosen & = 1 so 

that fi = S/ (Stt) = 0.848826... The domain s < ^ corresponds to i? > where the explicit form of is known (Eq. (|90|l 'l. 
^{s) is a smooth function with a very weak non-analyticity at s = ^ (essential singularity) -that can not be seen in a simple 
plot of $(s). 



2. Case R<0 (s > fi) 



The previous regime (density with support over ]0, L]) is only valid for i? > or equivalently s < fi. When R < 0, 



the effective potential is indeed minimal for x — xq 



> 0: the density pc is expected to have a finite support over 



[Li, L2] with Li > 0. Li and L2 are fixed by the constraints Pc{Li) = = Pc{L2). 

In this subsection, we find an expression for pc when i? < as a sum of elliptic integrals. We also derive the 
equations associated to the constraints pc{Li) = = pc{L2)- But we could in general neither compute explicitely 
the constraint J ^/xpc{x)dx = sVb nor find a closed form for the energy (and $), except for the asymptotic regimes 
s +00 and s For s p,'^ , we show that ^(s) has a very weak non-analyticity -an essential singularity- at 



p: 



$+(s)-$-(s)w-7r\/6(s-/i)e "^^(-^.) e'^C"^-!) 



s — > /i^ where ^{s) 



$ (s) for s < p 
$+(s) for s > p 



(97) 



The (normalized) solution pc, with support over [Li,L2], of the integral equation ([53]) is again given by Tricomi's 
theorem. We get 



with 



1 



TT^y/x — Liy/L2 — X 



-+-^^L,+L2-2x) + 



(t-y) VtTl 



2^1 -y) ^ + y 1 



e+y Vi+U (e+y)VT+e V2/(i+e)'i+e 



(98) 



(99) 
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where K and E are the complete elhptic integrals of the first and second kind respectively ; and 11 is the incomplete 
elliptic integral of the third kind: 







E{k)^ I J \ dt and K{k) ^ [ J- — dt (100) 
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FIG. 9; Density of states Pc{x) (density of charges) of the Coulomb gas associated to the computation of the pdf P{Gn = s-s/iV) 
of the center of mass, in the case s > /i = {R < 0), plotted for different values of the Lagrange multiplier R, or equivalently 

different values of s (and for 6 = 1). The effective potential seen by the charges is minimal for x = xo = > 0, thus the 
density has a finite support over [Li, L2] and is maximal around x = xq. 

When s > and s increases (i.e. the center of mass is larger than its mean and increases), L2 > 4, Li > and L2 and Li 
increase also: the charges form a bubble that gets further from the origin when R decreases (or s increases). 



We want to compute the pdf P{Gm = sv N). The basic variable is thus s. There are now four unknown parameters: 
R and D are two Lagrange multipliers and Li and L2 are the bounds of the density support. These parameters will 
be determined by enforcing the four constraints p{x)dx — 1, ^/xp{x)dx — s Vb, Pc{Li) = and Pc{L2) = 0- 

For a given R (Lagrange multiplier), the parameters Li and L2 are fixed by the constraints Pc{Li) = = pc{L2)'- 

r— RK{k) , , L2-ii , ii 

■\/L2 = where k = = 1 — — 

TT L2 L2 

and |J = -Kik) (^E{k) + - l) Kik)^ (102) 

We have already taken account of the constraint J pc{x)dx ~ 1 (normalization) by setting the constant Co that 
appears in Tricomi's theorem (equation (j85p ) to Cq ~ tt J Pc{x)dx — vr. 

The last constraint / \/xpc{x)dx — s^fb gives i? as a function of s. But the integral is in general difficult to 
calculate. And finally D is in principle given by the saddle point equation (see Eq. (I82|) ) at a special value of x, for 
example x = Li. But it is again difficult to compute in general. 

Therefore we couldn't compute exactly the saddle point energy. But, thanks to the above formulas, we could plot 
the density for different values of k (or equivalently, different values of R or of s). In this phase (i? < 0), as figure [9] 
shows, the Coulomb charges accumulate in a band near the minimum of the effective potential. They form a bubble 
that gets further from the origin when R decreases (or s increases). In this case, the interfaces are not bound to the 
substrate. 

We could also derive the asymptotics of <&(s) in this regime: s +00 and 5^/1+. 



Right tail of the pdf: limit s +00 (i? —00) 
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The limit i? — s- — cx) or equivalently s — > +00 corresponds to L2 — > +00 with k = — s- 0"*". In this limit, we 
have 



L2 « f -f-| + o(fc) 



i? w V2 

And finally, for fc = 1 - ^ 0+ with y = ^^^l^ fixed, < y < 1, we have: 



and <^ as fc = 1 ^ ^ 0+ (103) 

ii « i-f + f + o(fc) ^2 



a; — ii 



« — - y) A: + 0(fc2) with y=- ^ (104) 

The constraint J y/xpc{x)dx — s \fh gives s ~ f \/f + 0(1) as fc — + 0+, and the minimal energy diverges: 
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^JPc] « TT + 0(-) as fc ^ 0+ thus $(s) w 5^6 + 0(5) as s ^ +cx) (105) 
k 

which corresponds to a Gaussian tail: 

P(GAr = sTiV) cx e"^^'^' as s ^ +c» (106) 

3. Non-analyticity of the pdf: limit s — > fi'^ (R 0~) 

In this subsection, we analyse the limit s which corresponds to i? ^ 0^. 

Let us define for convenience the following parameters: 

i = ^ ^\ (^ as s ^ ^+) and X = -^-^ + ln2 {X +00 as s ^ ^+) (107) 
L2 — Li 4 

In the following, f is chosen to be the small expansion parameter (for s (X^). We will see that the expansion terms 
are of order 0{X'' = O (| Inf C«) with 6* > 0. As X''^ > X"^' £,^' (|ln^|''^^ > lln^l'^'^*') for < 61 < 61' and 
for every r] and 77', we can make an expansion in powers of ^ of the form X^eoo C6i(A^) ■f^, where the exact value of the 
coefficients cg{X) can be computed as functions of X without expanding them. We thus keep all the orders of the 
expansion in X (expansion in In^). 

We will show that the saddle point energy (and thus the pdf of the center of mass P{Gn = s^/N)) has a very 
weak (infinite-order) non-analyticity at s — ^ (mean of the center of mass). More precisely, we will show that the 

diff'erence of the energy on the right and left side of is of order O (^-^^ ~ O jt^) ~ O {\s ~ 
an essential singularity (it is much smaller than any power of \s — ^\). 

A singular limit for the saddle point density 

Using the equations (jl02p obtained by enforcing the constraint pc{Li) = = pc(^2), we can expand the Lagrange 
multiplier R and the bounds Li and L2 in terms of the small parameter ^ = -j^r^i to first order in f: 



e ^Vbis-^i : it is 



^X{X-1) V 16 [A(X - 1)]^/' / 



R » , ^ + e ^ T7T ]+0[^] and <^ (108) 



with f = T-^- and X = -ii^ + In 2. 



For s we have ^ ^ (AT +00) and we recover R ^ (with i? < 0), L2 ^ 4 and Li 0. These are 

the same limits as on the left side of fi: for s p~ ,we have R ^ and the density has a support over ]0, L] with L ^ 4. 

The saddle point density is given by Tricomi's theorem in equation (|98p . Using the constraint pc{Li) = 0, we get 



p,(x) = p{y) ^Aj-^ + + 4== ^iti^ y = (0 < y < 1) (109) 

- y) - y) L2-L1 
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where J(^, y) can be expressed as the principal value of an integral (see equation ([99|l V 

and where the coefficients A — —t^, B = , ^.f'/}^'^} and C — . ., can easily be expanded to first order in f. 

For y e]0, 1[ fixed and for ^ (s ^ we have 

J(e,y)«-2 + 2yT^argth(v/r^)-^ + 0(e) and ii « 0(0 as C (111) 
Therefore, to zeroth order in ^ the density shape (for Li < x < L2) is the same as for s < /i, it diverges for small x: 



But (for s > ij), the density /3c(2:) has a finite support over [Li,L2\ with Li > 0: it must vanish at x = Li. The 
constraint Pc{Li) = seems to be violated in Eq. (|112p . but it is not. As Li w 0{£), the part of the density associated 
to small X (close to Li) - and where the density must approach zero- does indeed not contribute to the zeroth order 
expansion of the density. The weight of the small range of values of x (around Li) where the density grows from zero 
to a very large value just becomes negligible when s — > /x"*". 

The limiting shape of the density for s ^ /i+ is thus singular. Therefore it is better not to expand J((^, y) and 
the density for fixed y (fixed x) and small ^ but to directly make an expansion of the energy, that involves integrals 
such that / dyJ{^, y)\/y + £, or J dyJ{^,y)hiy. Otherwise, as the limits y ^ and ^ ^ do not commute, the 
expansion of J(0 y) in terms of powers of ^ will generate increasing negative powers of y that will make integrals like 
/ dyJ(^, y) Iny diverge in zero. 

Expansion of the constraint J dxpc{x)^/x = s\/b for s ^ 

We must enforce the constraint dxpc{x)^/x — sVb that replaces the delta function S ^^^n J]^^^ ~ s\/lfj in 
the expresssion of the pdf of the center of mass P{Gn = s\/N): 

sVb = ( ' dxpc{x)^^{L2-Lif'^ ( dyp{y)^in (113) 

J Li JO 

From the expression of pc given in Eq. (|109p . we see that we need to expand for small ^ a double improper integral: 



m = dy^l^dM = r dy^M= V r dt (114) 
io yW^) Jo 'VvO^) Jo Vt+l t-y ^ ' 

As /(i^) is a double improper integral (with principal value), it is not easy to compute it or even expand it directly 
(for small Let us first make a simple transformation in order to get rid of the principal value; 



m = m - 0) +Uo{0 with /o(0 = dy f dt ^Jj—X r , ^ n (115) 

Jo Jo vy(i - y) VtT^ [VKy + + Vvit + f ) 

where J(f = 0) = —2 (it can be easily computed exactly) and where /o is a definite double integral, easier to expand. 
However, as we already noticed, the limit — *■ and the integration do not commute: the expansion can not be done 
inside the integral. Hence, the method of expansion must be a bit more subtle. Our method (see appendix-B for 
details) consists in splitting the initial integral /o(C) in ^ sum of integrals (some of them are easier to compute, the 
other ones are shown to be negligible). 

Finally (see appendix-B) we get the expansion of /(^ to first order in ^ (but to all orders in X, or In^: 

K -2 + ^\8X^ -AX -l\+0{^^ X^) as^^O (116) 

Hence the constraint J dxpc{x)^/x = s^/b is given by 

r- 2(4:X-3)VX /-16A:3 + 12A:2-2X + 1\ , ^ , , 

sVb K —— -TT7^ + ^{ ^^]+0(fX) asC^O (117 

3^ X- 13/2 1^ Sn{X-l)5 ^VX y V', y s \ ) 
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In particular, as expected, when ^ — > 0"'" {X ^ oo), s Vb tends to the mean value fj,Vb = 

Moreover, the formula above (Eq. pi7p l can be inverted to express X and ^ as functions of (s — /i). As /i \/b = 
we get: 

^ + l + 0((s-^)) and ^«e-(=-')^^ e^^^"^-!) (i + O ((s _ ^,))) as s ^ Ai+ (118) 
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Energy i?s[Pc] and scaling function (f>(s) = i?s[Pc] 2 
From equation (jSOp . we can compute the saddle point energy: 



1 /"^^ /2 Y) 

EsIPcI^t; dxpc{x)x- —sVb~ — (119) 



2 J^^ ' 2 2 

where the Lagrange multiplier D can be calculated by replacing x by Li in the saddle point equation for the density 
(equation (|82p ) and where J^^^ dx pc{x)x is not very difficult to expand for ^ ^ 0. Finally we get the expression of 
the energy for s ^ /i+ (^ ^ 0) , to first order in ^ and all orders in X = In 2 — ^ : 

r 1 . (X-l\ {iX^-AX + 2\ /-16X3 + 16X2-6X + 1\ 

[^^1 - ^H^) + I 2(X-1)^ ) + ^ ( 8X(X-1)3 j + (120) 

Using equation (|118p giving the expression of ^ as a funtion of (s — /i) in the limit s — s- /i+ , we will thus derive the 
behavior of the pdf of the center of mass P{Gn = s^/N) cx e^^ E^lp^] (^f^j. j^rge N) for s — > 

In order to show that the pdf of the center of mass has a non-analyticity at s = /i, we must compare the expansion 
of the saddle point energy on the left side and the right side of the mean. 

Zeroth order in ^: $(s) seems to be a smooth function 

Let us first consider the zeroth order in the expansion in terms of powers of ^ (on the right side of fi) . To zeroth 
order in f , the constraint f dxpc{x)\/x — sVb given in Eq. (|117p reduces to 



^3/2 j-1/2 

3n{X - 1)3/2 ' --y^' - 



.^„i<H^+o«,.#:+^+o(o ,121, 



Therefore, to all orders in In^ (or X -but to zeroth order in ^), we recover the same equation as Eq. ([57)1 . i.e. the 
same equation as on the left side of the mean, giving L {L2) as a function of s! The Lagrange multiplier R is also 
given, to zeroth order in ^ by the same function of L2 (L) as on the left side of the mean (see Eq. (|57)) ): 

+o(e)«i? + o(e) (122) 



VL^ 2 ^X{X-1) 

Finally, the energy, to zeroth order in ^ (but to all orders in X or In^) is given by the same expression as the energy 
on the left side of the mean: 

r 1+ , fX-l\ fSX'^ -4X + 2\ ^ 



2 In ^-^ + 1 + O , 
32 \ 2 ) \X 

« Es[pc]-^ol^^ (123) 
where L2 — ^2(5) is given by equation (|12ip . the same equation for s — > p.'^ to zeroth order in ^ as for s — > p~ . As 

-8 

^ fn Q7vis-^)s/b g4(in2-i) .^yjj^gj^ s /i+ (see equation (IllSp ). we get: 

$+(s) - $-(s) = E4pc]+ - Es[pc]- « O (^^^ « O (|s - Ml e"i=-'i^) (124) 

All the terms of the expansion of the energy (and thus ^(s) and the pdf of the center of mass) in powers of \s — p\ 
(or or i) are thus the same on the left and right side of the mean: $(s) is a smooth function, it is infinitely 
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difFerentiable even at s = /i -in particular the quadratic approximation of ^{s) in Eq. ((91]) is valid on both left and 
right side of its minimum (s = /i). However, we will show that the expansion to first order in (by keeping all the 
powers of X) gives a very weak non-analyticity of the energy (and thus 4'(s)). 

First order in ^: non-analyticity of $(s) 

Using equation (I120|) and the remarks we made about the zeroth order expansion in ^, we get the difference between 
the expansion of the energy on the right and left side of /i: 

Using the expression of ^ and X = In 2 — as function of s for s /i+ given in Eq. (jllSp , we finally get 

$+(s)-$-(s) ^ Esipc]-^ - Es[pcr w -7r\/5(s-^)e~'vi;t-M) g^O^^-i) as s ^ /i+ (126) 

This is an essential singularity. We have shown that the pdf of the center of mass P{Gm = s VN) oc E^lp^] ^ 
very weak non-analyticity at s = fi: the energy (or equivalently *&(s)) has an infinite-order non-analyticity, of order 

O (^\s - p.\ e~^^^^^^^ 



IV. CONCLUSION 



In summary, we have studied a simple model of TV nonintersecting fluctuating interfaces at thermal equilibrium 
and in presence of a wall that induces an external confining potential of the form V{h) = ^ — h "^"^^^-^ . Our study 
extablishes a deep connection between the statistics of heights of the interfaces in the limit of a large system {L oo) 
and the eigenvalues of the Wishart random matrix, thus providing a nice and simple physical realization of the Wishart 
ensemble. More precisely, we have proved that the joint probability distribution of the interface heights hi in the limit 
of a large system can be mapped to the distribution of the eigenvalues Xi of a Wishart matrix under the change of 
variables 6 /if = , with arbitrary parameter value M — iV of the Wishart ensemble that is fixed by the parameter a 
of the inverse square external potential. 

We have also shown how to exploit the relation between interfaces and eigenvalues of the Wishart matrix to derive 
asymptotically exact results for the height statistics in the interface model. In particular, we have seen that the 
nonintersecting constraint, the only interaction between interfaces in our model, drastically changes the behavior 
of interfaces: they become strongly correlated. Despite the presence of strong correlations that make the problem 
difhcult to analyse, we were able to compute a number of asymptotic (large N) results exactly. These include the 
computation of the average density of states, the distribution of maximal and minimal heights and the distribution 
of the center of mass of the interfaces. In the last case, we have shown that the distribution has an extraordinarily 
weak singularity near its peak (an essential singularity) and this non-analytical behavior was shown to be a direct 
consequence of a phase transition in the associated Coulomb gas problem. 

Finally, we expect that the appearence of the Wishart random matrix in a physically realizable example as shown 
in this paper will be useful in other contexts. In addition, the Coulomb gas technique used here seems to be a very 
nice way to derive exact asymptotic results in this class of interacting many body systems where exact analytical 
results are hard to come by. It would be interesting to use the analogy with a Coulomb gas in other physical problems 
related to Wishart matrices, for example to compute the distribution of entropy of a bipartite quantum system (see 

EES)- 

Acknowledgements: It is a pleasure to thank A. Comtet for many useful discussions. 

APPENDIX A: COMPUTATION OF THE MOMENTS OF THE MINIMAL HEIGHT 

([50)1 gives an exact expression for the pdf of the minimal height (lowest interface): 

P ( Vin ^t,N) = 2b^ e-^^*' (-& t') = e-^^*' 7V(7V + 1) iFi (1 - iV, 3, -6^^) (Al) 

(see ([14]) for the relation between Laguerre polynomials and hypergeometric functions) 
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Therefore we can compute explicitely the moments of the minimal height: 

r°° f°° 2 

(^min) = / dtt'' P{K-,^=t,N)^b'N{N + 1) / dtt^+'^e-^^' iFi{l- N,3,-bf) (A2) 
Jo Jo 

lj-k/2-2 poo 



= b^N{N+l)' 



du g-JVn ^^^(1 _ JV, 3, ~U) 



(A3) 



with bt^ — u. The integral above can be computed: duu'^ ^ e iFi(a, b, — u) = c r(d) 2^1 (a, d;b;—l/c) 
Therefore 



{h 



K ) 

mm/ 



r(fc/2 + 2) (iv + i) A. i./9j_9'? i//\n 

2;,fc/2 Nk/2+1 2i^i(l-iV,fc/2 + 2;3;-l/7V) 



For example, for fc = 1, we find: 



r(5/2) (jV + 1) 

2fel/2 7V3/2 



2i^i(l-iV,5/2;3;-l/7V) 



For large N, as ^^^^(l - A^, 5/2; 3; -1/A) = E„ ^^{gKS^+g^' (i-A-)(2-jV).,„-A.) ^ ^^^^ 



lim 2Fi(l- A,5/2;3;-l/Ar) = V 

N^oo ^ — ^ 



(5/2)(7/2)...(3/2 + n) 
(3)(4)...(n + 2)n! 



Hence, for large A^: 



with 



= iFi(5/2;3;l) = ^/o(l/2) 



Cl 



= ^¥ = #/o(l/2) « 1.5538 



(A4) 



(A5) 



(A6) 



(A7) 



APPENDIX B: NON-ANALYTICITY OF THE PDF OF THE CENTER OF MASS: EXPANSION OF 1(0 



Let us expand for f the integral !{£,) given in Eq. ([TT4]): 

,,,, 



(Bl) 



As /((^) is a double improper integral (with principal value), it is not easy to compute it or even expand it directly 
(for small 1^). Let us first make a simple transformation in order to get rid of the principal value: 



no 



\/y{i-y) Jo Vt + l t-y 



-2 + e i'dy fdt 



^2/(1 - y) ^-y VV 

1 



/O ^0 

= -2 + Uo{0 
The value of /(^ — 0) can indeed be computed exactly: 



^2/(1 - y) Vt + ^ V^y + I) + VW+I) 



V dt-^ ' ^ 



= 0) = / dy- , , , , 

/o v^-y Jo t-y /o V 1 - y 



-2 + 2Vl-yargth (^1 



(B2) 



(B3) 
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/o is a definite double integral, easier to expand: 

/o(0= fdy fdt^^. 
Jo Jo vy(l - y 



(B4) 



We thus need to expand a definite double integral /o(^) (we have got rid of the principal value). 

However, as wc already noticed, the limit ^ ^ and the integration do not commute: the expansion can not be 
done inside the integral. 

Let us thus consider separately the integration over ]0, ^] and 1[ (for the variable y): 



/o(o = /i(e) + ./2(o 

where /i and /2 are definite double integrals (no principal value): 



h{0 
h{0 



dy 



Jo Vy(i - y) .Jt + l yt{y + + Vy(i + 



dy 



Vvi^ - y) ^/t + l ^t{y + + Vy(* + 



(B5) 

(B6) 
(B7) 



The method of expansion will be the following. We will write /i (resp. /2) as a sum of two integrals: the first will 
be chosen of the same order of /i (resp. J2) for small ^, but easier to compute (by separation of the variables t and y 
for example); the second will be much smaller than the first and than /i (resp. f^)- Then, following the same scheme, 
each of these two integrals can again be split into two pieces, until we get the full expansion to first order in ^ (the 
last integral will be shown to be much smaller than the other and will be neglected). 

Let us first make the change of variable y = in /i(^): 



/i(0 



d'li 



dt 



1 



v^yr^i^ViT^ (y^i(rT^+0i(rro) 



Then we have 

/i(C) = ./f(0 + /i'(0 

with f1 a product of two integrals (separation of variables) 



(B8) 



(B9) 



^/u {y/u + Vi + m) VI - 

du 



^/u{^/u+ VT+u) 



+ 0(0 



a" 



dt 



(-1 + \/2 + argsh 1 + 0(0) (- In^ + 41n2 - 2 + O (^ In^) 

In^ (1 - 72 - argsh 1) + (41n2 - 2)(-l + V2 + argsh 1) + O (^ In^ 



(BIO) 



and /i = /i — fi (expected to be much smaller than /i) is given by: 



= [ du [ 
Jo Jo 



du 



^U{y/U+VT+U) 



f 

Jo 



(V^(1 + m) + ^^(t + O) ViiVu + y/l + u) 
dt- 



VWTo iVt + Vt+^) (^VW+uj + ^/WTo) 



fm + /i'(o 



(Bll) 



with 



dz- 



w ln2 (^3\/2 - 4 - argsh + In^ 
where we have made the change of variables t = ^z; and 



Jo VT^^{^/^ + Vmi)^ Jo + (Vi + ViTI)' + + + «)) 



/■I /-CX) 

du 
Jo Jo 



+0(0 

« argshl(ln2-3) + V2(l-31n2)- l + 61n2 + O(0 
Thus we have, for ^ 

/i(0 « In^ (l - V2-argshl) + (^1 - + argsh 1 (41n2 - 5) + ln2 (4\/2 - 2)) +O(^ln0 

The same method of expansion aplied to /2 gives 

/2(0=/3(0+/4(0 

with 



M)!+ln^ (^-31n2) +21n2 (21n2-0 +O(^(ln0^) 



and 



dt. 



/4(0 = / ■ — , / , . 

Ji yVT^{Vy + vy + ^) Jo .jt + l {Vi + Vt + I) [Vt{y + + Vvit + O) 

w 1^-^ -ln2 + \/2 + argsh 1^ In ^ 

+ ((5 - 4 In 2) argsh 1 + 4(ln 2f + In 2 (1 - 4\/2) -2 + ^2^+ 0(C(ln ^f) 
(for the expansion of /4, the same method of splitting has again been applied) Hence 

/2(0 = /3(0 + m) « + In? (\/2 + argsh 1 - 41n2) 

+ ((5 - 4 In 2) argsh 1 + 8(ln2)2 + In 2 (-2 - 4^2) - 2 + V2j + 0{^{\n^f) 
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and, as /(C) = -2 + ^ + /2(e)), we get (with X = \n2- 



m 



(InC)^ 



InC (1 - 4 In 2) + (8(ln2)2 - 4 In 2 - l) + 0(^2 (in^)2) 



-2 + ^ 8x^-4^-1 +o{s;x 



(B19) 
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